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are easily satisfied. However, the model is well constrained by the indirect detections of 
VDM from BBN, CMB, AMS-02, and diffuse y/X-rays. Consequently, only when the dark 
Higgs boson mass is at most of O(keV) does there exist a parameter region which leads 
to a right amount of VDM relic abundance and an appropriate VDM self-scattering while 
satisfying ah other constraints simultaneously. 
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1 Introduction 

In spite of increasing astrophysical and cosmological evidence for the existence of the dark 
matter (DM) [1, 2], the nature of DM remains a mystery. According to the dominant 
paradigm DM consists of collisionless, cold particles that successfully explain the large 
scale structures in our Universe. However, collisionless cold DM predictions obtained by 
N-body simulations face some difficulties known e.g. as the cusp-vs-core problem [3-6] or 
the too-big-to-fail problem [7-9] when confronted with precise observations at the dwarf 
scale. However, it has been shown that the presence of sizable DM self-interactions with 
(Jdm/^dm = 0.1 ~ 10 cm 1 2 /g has the potential to alleviate such a tension [10-17], even 
though the DM self-interactions are constrained to be <tdm < 1 cm 2 /g by measurements 
at the cluster scale [18-23]. 

Such large DM self-scatterings naturally arise if there is a light particle mediating 
the DM interaction and the corresponding cross-section is enhanced by non-perturbative 
effects [24-33]. One immediate consequence of this light mediator scenario is that the 
DM self-interaction cross section is velocity dependent [14-16, 19-22], which allows for 
the signals at the dwarf scale to evade the constraints from the galaxy clusters. A simple 
way to realize this scenario is to introduce a model, where DM is generated via the dark 
freeze-out mechanism in which it predominantly annihilates into a pair of light mediators. 
Nevertheless, it has recently been shown in refs. [34, 35] that this secluded DM model [36] 
is severely constrained by the DM indirect detection. A way to avoid these problems is 
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to consider a DM production mechanisms different from the conventional freeze-out. One 
possibility is the freeze-in mechanism [37, 38] (see i.g. ref. [39] for a recent review and the 
complete references therein). It is found in refs. [37-41] that the final DM relic density is 
determined exclusively by the main DM production channels at the freeze-in temperature 
and it is not sensitive to many details of DM evolution at higher temperature, which 
guarantees the predictability of this mechanism. 

The freeze-in as a production mechanism for self-interacting dark matter was analyzed 
in [42-47]. Notably, the case of light-mediator was discussed in ref. [44] within the model of 
Hidden Vector DM with dark SU(2) gauge symmetry [48], where it has been found that the 
scenario with keV mediator agrees with experimental constraints. It has been also noticed 
that if decays of the mediator into e + e - are allowed, its significant abundance and large 
lifetime cannot satisfy bounds from Big Bang Nucleosynthesis (BBN), so that this region 
of the parameters is excluded. 

In this work, we study an abelian version of vector dark matter (VDM) models [48- 
58] in which the VDM particle with mass of 0(GeV ~ TeV) couples to the SM sector 
only through the Higgs portal. We take into account recent bounds from BBN, CMB and 
discuss possibility of constraining the model with FERMI-LAT, AMS-02, diffuse y/X-Ray 
and direct detection LUX data. In the case of indirect constraints on DM annihilation, 
we include the effect of Sommerfeld enhancement. We also take into account consequences 
of electro-weak phase transition in calculation of DM production. The dark Higgs boson 
of the VDM model is assumed to be so light that it can induce large self-interactions to 
solve the small-scale structure problems. We focus on the scenario in which the VDM is 
produced by the freeze-in mechanism. The main question that we address is if there exist 
a region in the parameter space that can generate the right VDM relic abundance and 
appropriate DM self-scatterings while still satisfying all the direct and indirect detection 
constraints. After scanning over the parameter space we conclude that if the mediator /12 
is too light to decay into e + e - , then indeed all the constraints can be satisfied together 
with correct relic abundance and appropriate DM self-scatterings. The necessary mediator 
mass is of the order of O(keV). Our results agree with those found in [44]. 

The paper is organized as follows. In section 2, the VDM model is presented. The 
VDM production via freeze-in is discussed in section 3, with a special attention to the effects 
of the electroweak (EW) phase transition. Then we discuss the VDM self-interactions in 
section 4. Section 5 and section 6 are devoted to constraints from DM direct and indirect 
detection experiments. The numerical results are presented in section 7. Finally, we give 
a brief summary in section 8. Some useful formulae are collected in appendix A. 

2 The model 

Following refs. [53, 54], we introduce a dark U(l)x gauge symmetry and a complex scalar 
S which is neutral under SM gauge group but has unit charge under this U(l)x symmetry. 
We further assume an additional Z 2 symmetry, under which the gauge boson X M and S 
transform as follows: 

X^-X^ ( 2 . 1 ) 
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which is just the charge conjugate symmetry in the dark sector. It forbids the kinetic 
mixing between the SM U(l)y gauge boson B M and , ensuring stability of X /JL . 

Therefore, the relevant dark sector Lagrangian is given by 

Cd = + (D^D»S+/4\S\ 2 - As|5| 4 - k\S\ 2 \H\ 2 , (2.2) 


where H is the usual SM Higgs SU(2)l doublet, and the covariant derivative of S is defined 
as D^S = (d^ + igxX^S with gx being the corresponding dark gauge coupling constant. 
Note that the quartic portal interaction, k\S\ 2 \H\ 2 1 is the only connection between the 
dark sector and the SM, so in the limit hi —)► 0 the two sectors decouple. Also, the mass 
term of S has the negative sign compared with the usual scalar field, so that it can induce 
the spontaneous symmetry breaking (SSB) of the gauge U(l)x- By minimizing the scalar 
potential of the model, we can obtain the vacuum expectation values of the usual SM Higgs 
doublet ( H ) = (0, vh/V%) T and the dark Higgs (S) = vs/V% as follows: 


_ 4As/i| 2^/ig 2 _ 4A hI^s 2hifi 2 H 

~ 4A h A s - ’ VS ~ 4A^A s - 


(2.3) 


Note that ( S) can be always assumed real without compromising any generality, therefore 
the discrete symmetry (2.1) remains unbroken as needed for the stability of X^. 

After the SSB happens, the dark gauge boson obtains its mass mj = gxvs y ia the 
dark Higgs kinetic term, and both scalar fields can be written as 


H = 


h+ ) . 

[v H + (f>H + IVh)/) 


S = -^( v s + <j>s + ivs) ■ 


(2.4) 


By expanding the scalar potential up to the second order, the mass squared matrix M 2 of 
the two physical scalars {4 > Hi4 > s) T is given by 


M 2 = ( 2:KhV h KVhVs ) 

\ kv h v s 2 A 5 u| J 

With the following orthogonal transformation of scalars, 


4>h 

<Ps 


( c e -sb 
yse ce 



(2.5) 

( 2 . 6 ) 


we can define the mass eigenstates ( h \, h 2 ) T with their masses (//</,,, ), where 6 is the 

mixing angle with sq = sin 0 and eg = cos 0. As a result, we have the following relations: 

m 2 hl -ml 2 )s 2 e m 2 hi c 2 e + m 2 h2 s 2 d m 2 h <%-+ m 2 hi s 2 

s = — 2^r s —' Xh = — 24 —• As = — 24 —- (2 - 7) 

In the freeze-in mechanism, the dark sector composed of X and /12 never thermalizes 
with the visible SM sector, so that the portal interactions hi or sq should be very tiny. 
As is evident from eq. (2.6), the h\ boson is mostly SM-Higgs-like, while is almost 
the dark Higgs <j>s. We have found that the most convenient choice of input parameters 
which specify the models is (rax, ra^ 2 , ft, gx), together with the already known parameters 
vh = 246 GeV and ra/ ll = 125 GeV. 
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3 Vector dark matter relic density via freeze-in 


Within the freeze-in mechanism, the standard assumption is that the initial abundances of 
the VDM and the dark Higgs after reheating are assumed to be negligibly small, which 
is possibly a result of the reheating process itself or another mechanism. Furthermore, the 
Higgs portal coupling should be very tiny so that the dark sector can neither thermalize 
itself nor be in equilibrium with the SM sector. When the VDM mass is smaller than the 
EW phase transition temperature Tew — 160 GeV [59, 60] its abundance is mainly con¬ 
trolled by various SM particle annihilations and/or decays that contribute to the collision 
term of the following Boltzmann equation 



(3.1) 


/ 


where Yx — ux/s is the DM yield defined as a ratio of DM number density nx and the 


entropy density in the visible sector s. The parameter x = mx/T describes the SM sector 
temperature T, H is the Hubble parameter, and 7 ^ = (av)i(nl q ) 1 2 denotes the so-called re¬ 
action density [41] for the SM particles annihilation into VDMs (for yj we sum over all SM 
fermions). The last term represents the SM-Higgs-like h\ decays to a VDM pair when this 
channel is kinematically allowed. Since here mh 2 <C mj, no corresponding decay term for 
/12 appears. In this project, the model is implemented within LanHEP [61, 62] and calcula¬ 
tions of the cross sections and decay rates are performed adopting CalcHEP [63]. Definitions 
of reaction densities, obtained cross sections and decay rates are collected in appendix A. 

It is interesting to note that all of the reaction densities are proportional to the square 
of the Higgs portal coupling ft 2 with no dependance on gx} which explains why we have 
decided to use ft instead of sin 6 as a parameter. Also, due to the assumed mass hierarchy 
rrih 2 <C mj, the value of m^ 2 influences the resulting DM abundance very weakly. Hence, 
the prediction for VDM relic abundance depends mainly on two parameters mx and ft. 
Figure 1 shows typical examples of evolution of reaction densities for various SM channels as 
functions of the SM sector temperature T. The left and right panels represent the case with 
mx larger or smaller than rn/^/2. It is evident that in the first case, only the annihilations 
of the SM particles are involved. On the other hand in the second scenario the decay of 
the visible Higgs h\ is open and overwhelms other annihilation modes near the freeze-in 
temperature Tpi . 2 Since the freeze-in mechanism is IR dominated [38, 41], the VDM relic 
density is dictated by the h\ —* XX decay rate. We present the resulting evolution of the 
VDM yields Yx in figure 2, which illustrate typical features of the freeze-in mechanism. 

However, when the VDM mass is much larger than the EW phase transition temper¬ 
ature Tew 5 the VDM abundance stops increasing before the EW phase transition. In this 
case, the SM gauge symmetry SU(2 )e x U(l)y is not broken, so that only the tree-level 
diagram shown in figure 3 can generate VDM particles. Hence, the Boltzmann equation 

1 This property could be easily seen adopting the first relation of (2.7). 

2 Defined as the temperature at which the VDM yield from freeze-in production reaches its maximal 
value and stabilizes. 
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m x = 100 GeV 


m x = 50 GeV 
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Figure 1 . Evolution of the reaction densities for various SM channels as functions of the SM 
sector temperature T for mx = 100 GeV and k = 2.09 x 10 -11 (left panel) and mx = 50 GeV and 
k = 2.40 x 10 12 (right panel). The chosen values of k result in the observed DM relic abundance. 


A 7? x =100 GeV, /c=2.09*10 -11 m x =50 GeV, /c=2.40*10“ 12 




Figure 2. Examples of evolutions of the VDM yield Yx as the functions of x = mx/T , in which 
the model parameters for both panels are chosen the same as those in figure 1. 



Figure 3. Feynman diagram to generate the VDM X via the SM Higgs doublet H annihilations. 


can be simplified to: 

X ^ S ~dx~ = ^ HS ' ( 3 - 2 ) 

By comparing the solutions to the Boltzmann equations in eqs. (3.1) and (3.2) with 
the observed DM relic density fly/i 2 = 0.11, we can obtain the value of k as the function 
of the VDM mass mx in figure 4. It is interesting to see the change of K-mx scaling in 
this plot, which reflects the transitions of the dominant VDM freeze-in channels. When 
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Figure 4. The value of k as a function of the VDM mass mj, that gives the observed relic density 
via the freeze-in mechanism. 



the VDM mass is larger than m^/ 2, as mentioned before, only the annihilation modes 
contribute, no matter whether the EW gauge symmetry is broken or not. In this case, 
according to the argument in refs. [38, 40], the yield could be estimated as 


Yx ~ a(T F i)Mp\T F i ~ k 2 —^ ~ , (3.3) 

T f i m x 

where the first relation follows from the dimensional argument with M F \ being the Planck 
mass. cr(Tpi) is the total cross section of the SM particle annihilation at the freeze-in 
temperature Tpi, which is simplified to be a ~ k 2 /T^. We have also used the relation 
T F i ~ mx, which can be understood as follows. When m x > Tew? as it has been mentioned 
above only the channel HH t XX contributes to VDM generation. It becomes ineffective 
as the temperature drops below m x , since then the SM Higgs doublets do not have enough 
kinetic energies. On the other hand, for the case with m x < Tew? the VDM freeze-in 
process is dominated by the annihilations of particles which are lighter than the VDM. 
Similarly, when the SM plasma temperature decreases below m x , the VDM yield ceases 
to grow any more due to the fact that these channels are no longer kinematically allowed. 
Concluding, the freeze-in temperature is expected to be around the VDM mass, Tpi ^ m x , 
in the present scenario. Then it is easy to derive from eq. (3.3) that the predicted VDM 
relic density fl x h 2 oc Yxmx should only depend on k, whereas the dependence on m x are 
cancelled out, which is manifested as a flat line in figure 4. However, if the VDM is lighter 
than a half of the visible Higgs mass, the decay channel hi —>► XX dominates, so that the 
VDM yield should be 




h^XX 


M F \ 

T 2 
1 FI 


K 2 m h 


T 2 
1 FI 




,Mpi 

m h i 




(3.4) 


JHEP01(2018)020 








where the decay rate should be Thi-+xx ~ m h 15 and the freeze-in temperature in this 
case is Tpi ^ at which the density of the visible Higgs h\ is greatly suppressed by its 
Boltzmann factor. Hence, the VDM relic density is ftxh 2 oc K 2 mx, which results in the 
scaling of ft oc m x 7 in figure 4. Finally, note that the small but abrupt rise of ft at the 
mx — 160 GeV represents the EW phase transition effect due to the sudden change of the 
main VDM production channels. 

In order for the freeze-in mechanism to work, it is required that the dark sector neither 
thermalize by itself nor with the SM sector. It is easy to check that the portal coupling 
ft implied by the VDM relic density is so tiny that it is impossible for the dark sector to 
equilibrate with the visible one. However, the non-thermalization of the dark sector by 
itself is not guaranteed. When the number densities of the VDM and accumulated via 
freeze-in become large enough, it is probable that the dark sector process XX —» / 12/12 
would be cosmologically efficient, which would soon change the number densities of VDM 
and /12 to form a dark plasma with a common (and in general different from the SM) 
temperature. Therefore, one should ensure that thermalization in the dark sector cannot 
take place and the appropriate condition can be coded by the following inequality [41, 44]: 

(a(XX h 2 h 2 )v)nx < if, (3.5) 

where (a(XX ^ 2 ^ 2 )^)? ^ 1 , and H represent the thermally averaged cross section 
for VDM annihilations into /12 pairs, the number density of VDM, and Hubble param¬ 
eter, respectively, all of which are evaluated at the freeze-in temperature Tpi. Note that 
(a(XX h 2 h 2 )v) is proportional to the dark gauge coupling so that it should not 
be suppressed in the parameter space where the DM has large self-interactions. Thus, the 
condition in eq. (3.5) is not easy to be satisfied in the present scenario and therefore it 
constraints substantially the freeze-in parameter space as shown below. 


4 Vector dark matter self-interactions via a light mediator 


It is well known that the cosmological small-scale structure problems, such as the ‘cusp 
vs. core’ and the ‘too-big-to-fail’ problems could be ameliorated if DM self-interaction was 
sufficiently strong at the dwarf galaxy scale [10-17], the required value of the cross-section is 


^ „ cm 2 <jt . ~ cm 2 

0.1 - < — < 10 -, 

g m x g 


(4.1) 


where ctt = f dfi( 1 — cos 6)da/dft is the so-called momentum transfer cross section between 
DM particles. However, DM self-scattering cross-section as large as ax/mx — 10 cm 2 /g 
is not allowed by observations at the cluster scale with the typical constraint ax/mx < 
1 cm 2 /g [18-22]. 

A possible strategy that may generate large DM self-interaction is to introduce a 
mediator which is much lighter than the DM particles. In the VDM model, the elastic 
DM scattering is mediated by an exchange of the two Higgs scalars, h\ and / 12 . In the 
limit of small mixing, the /ii-mediated contribution is negligible due to sin a and large hi 
mass suppression. In contrast, XX /12 coupling is not suppressed by small mixing and, 
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in addition, it is much lighter than the VDM particle, therefore can act as a light 
mediator which might be capable to amplify the self-interaction. When axmx mh 2 
with ax = g\/{ 47 r) the fine-structure constant in the dark sector, the perturbative Born 
approximation is applicable in which the dominant t-channel /^-exchange to the transfer 
cross section as follows [29]: 


^Born 

(Jji 


87 ra x 

TTl^V 4 


In 



2 2 
m x v 


m 


h 2 


2 2 
m x v 

2 I 2 2 J 

mji 2 + m x v z 


(4.2) 


where v is the relative velocity in the VDM two-body system. Nevertheless, beyond the 
Born range, /12 is much lighter than axmx so that the nonperturbative effects would 
become important and give rise to the following attractive Yukawa potential: 


V{r) = 


ax& 


-m h2 r 


r 


(4.3) 


Note that due to such nonperturbative corrections, the DM self-interactions have the non¬ 
trivial dependence on the VDM velocity. When the range of the potential characterized by 
l/rrih 2 is much larger than the VDM de Broglie wavelength 1 /(mxv), i.e., mxv m^ 2 , 
this part of parameter space is well known as the classical regime, for which analytic fitting 
formulas for <jt [28, 29, 31, 64] are available in literature. In our numerical calculations, 
we adopt the more recent improved analytic expressions provided in ref. [31]. On the other 
hand, if mxv < m^ 2 , the VDM self-scatterings can be enhanced by several orders of mag¬ 
nitudes due to the formation of the quasi-bound states. This region of parameter space 
is usually denoted by the resonant regime. In this work, we obtain gt in this regime by 
closely following ref. [29] to solve the non-relativistic Schrodinger equation with the poten¬ 
tial in (4.3). Moreover, it has been found [28, 29, 31, 64] that, with the presence of the non¬ 
perturbative effects, the VDM transfer cross section gt is enhanced more significantly as the 
relative DM velocity becomes small. Such a velocity dependence of VDM self-scatterings is 
very appealing, since it helps the VDM model to solve the small-scale structure problems 
for the dwarf galaxy scale with a typical velocity v ~ 10 km/s while evading the strong 
constraints from the galaxy clusters with v ~ 1000 km/s. More recently, a more careful 
analysis of DM self-interactions from a light (pseudo-)scalar has been presented in ref. [35], 
where a more appropriate definition of the momentum transfer cross section gt is given 
and the possible correction from the ^/-channel light mediator exchange is investigated. 
However, it is seen in ref. [35] that such corrections lead to very small modifications in final 
results so that we neglect them and follow the conventional formula from refs. [28, 29]. 


5 Direct detection of the vector dark matter 

It is usually claimed that the DM direct-detection experiments do not provide relevant 
constraints for models in which the DM particles are mainly produced by the freeze-in 
mechanism since the DM nuclear recoil cross sections are suppressed by tiny portal cou¬ 
plings. However, in the present scenario, the spin-independent (SI) VDM-nucleon (XN) 
scatterings are mediated by the two neutral Higgs bosons /ip 2 , and thus it is possible that 
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the cross-section is greatly enhanced by the small mass of the light mediator /i 2 - This 
feature is clearly reflected by the corresponding formula for the differential cross section of 
the XN scatterings with respect to the momentum transfer squared g 2 , 


dcrxN _ &xn 2\ 

dq 2 “ 4 / 4 ^ 


(5.1) 


where v is the VDM velocity in the lab frame, hxn = tuxtun/{ jnx + tux) is the reduced 
mass of the XN system, and 


&XN 


K2 fN m X m Nt i XN 

Krn^ml^m^+ALixNV 2 ) 


(5.2) 


is the total cross section for the XN scattering with the effective nucleon coupling fx ~ 
0.3 [65-67]. Compared with the usual definition of the SI independent DM-nucleon cross 
section in the literature, eq. (5.1) has an additional form factor G(q 2 ) defined as 


G{q 2 ) 


+ 4 ^JV^ 2 ) 

(q 2 +m 2 h2 ) 2 


(5.3) 


which encodes the effects of the light mediator /i 2 - It is clear that, for the heavy mediator 
case with m^ 2 q 2 ~ \(j? XN v 2 , the factor G(q 2 ) will be reduced to 1, i.e., we will recover 
the conventional XN contact interaction, and the usual experimental constraints can be 
applied. But when m^ 2 <C g 2 , the XN differential cross section in eq. (5.1) will have 
extra q 2 dependence characterized by the G(g 2 ), thus modifying the corresponding nuclear 
recoil spectrum and, in turn, the final fitting results. Therefore, we need to re-analyze the 
experimental constraints in the latter case. 

The strongest constraints on the direct detection of the VDM come from the LUX [68], 
PandaX-II [69] and XENON1T [70]. In the present work, we use the LUX 2016 dataset 
as an illustration of the SI direct detection limits to the VDM model since PandaX-II and 
XENON1T datasets would give the similar results. Due to the modification of the DM 
nuclear recoil spectrum caused by the light mediator / 12 , we follow the simplified analysis 
methods presented in refs. [71, 72] to give the LUX 90% C.L. upper bounds on the VDM 
nuclear scattering cross section <jxn and on the Higgs portal coupling ft, with the final 
numerical results as shown in figure 5. 

It is seen from figure 5 that the LUX upper bound on the VDM nuclear scattering 
cross section increases with the decrease of the mediator h<i mass. The lowest curve with 
mh 2 = 200 MeV corresponds to the point-like contact VDM nuclear interaction, and agrees 
with the upper limit in ref. [68], since such a mass of /12 is already much larger than the 
typical momentum transfer scale q ~ 10 MeV. However, when transformed into constraints 
on the Higgs portal coupling ft on the right panel of figure 5, the LUX upper bound is found 
to behave oppositely, that is, it becomes stronger with the smaller m\ l2 . Furthermore, 
when rrih 2 < 1 MeV, the LUX bound is shown to saturate a limiting curve, which can 
be understood that the /12 mass is cancelled out in the final expression in eq. (5.1) in 
this parameter region. However, even though it is remarkable that the LUX upper limit 
of ft reaches the order of lO -10 for large VDM masses, it is not able to give a meaningful 
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LUX bounds 


Parameter space for correct relic density and LUX bound 



Figure 5. LUX upper bounds on the total SI VDM nuclear recoil cross section cjxn (left panel) 
and Higgs portal coupling k (right panel) for different light mediator I 12 masses. The solid black 
curve in the right panel corresponds to parameters which reproduce the measured value of the DM 
relic density. 

constraints on the freeze-in region of our model. Thus, in the following, we will not consider 
the direct detection constraints any more. 

6 Indirect detection constraints on vector dark matters 

Phenomenology of indirect detection of VDM crucially depends on properties of the me¬ 
diator h, 2 , such as its mass m^ 2 , lifetime r^ 2 and dominant decay channels. Since we 
are interested in the light h<i which could give rise to the large enhancement of VDM 
self-interactions, we will limit ourself to m^ 2 < 100 MeV. Thus, the parameter space is 
naturally divided into two regions: (i) m^ 2 > 2 m e and (ii) rrih 2 < 2 m e , where m e is the 
electron mass. In the former region, the dominant /i 2 decay channel is e + e - pairs, while 
only the diphoton mode is kinematically available in the latter case. Consequently, the light 
mediator lifetime r^ 2 is different in these two regions. Specifically, 10 4 s < r^ 2 < 10 12 s 
in region (i) while r^ 2 > 10 12 s in region (ii), which is illustrated in figure 6 for a typical 
VDM mass mj = 100 GeV and a Higgs portal coupling k — 2.09 x 10 -11 consistently with 
the DM relic density (see figure 4). Analyzing constraints from DM indirect searches, we 
will consider these two regions separately. 

Region (i): 2 m e < m^ 2 < 100 MeV. Since dominantly decays into e + e - pairs 
resulting in its lifetime of 10 4 s < r^ 2 < 10 12 s, the relevant indirect detection constraints 
involve the following experiments: 

• Big Bang Nucleosynthesis (BBN): due to its long lifetime r/ l2 > 10 4 s, the dark 
Higgs boson h?, was still present in the early Universe at the epoch of BBN. Thus 
from the viewpoint of BBN, /12 plays the role of an extra, decaying component of DM. 
Such a late decay of would produce e + e - pairs with sufficient energy that would 
spoil the predictions of abundances of various elements [73-76]. We adopt the most 
recent results from ref. [76] where the authors also studied the BBN effects triggered 
by decays of dark Higgs bosons produced by the freeze-in mechanism. It is seen from 
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m x =100 GeV, k= 2.09 x 10 _11 



Figure 6. Contour plot for the ft, 2 lifetime in the m^-OLx plane for VDM with mx = 100 GeV. 
The Higgs portal coupling k = 2.09 x 10 -11 is determined by the VDM relic density through the 
freeze-in mechanism. The numbers on the line represent the ft 2 lifetime in units of seconds. 

figure 8 in ref. [76] that the most stringent constraint to the model is set by sq < 
5x 10 -12 for 1 MeV < m^ 2 < 100 MeV where 6 is the mixing angle defined in eq. (2.6). 
On the other hand, when m^ 2 < 1 MeV, there is no constraints to the decaying /12 at 
all. Note that the result in ref. [76] was obtained in the limit of k —>► 0 and vs —» 00 
while keeping 6 fixed, so the 2 —>► 2 processes involving top quarks predominate the 
production via freeze-in. However, in our scenario, the Higgs portal coupling does not 
approach zero. The most important contribution to ft 2 density arises from the SM-like 
Higgs decay h\ —>► ft2^2, which is more efficient than the top quark annihilations and 
top-gluon inelastic scatterings. Therefore, we expect that ft ,2 is more abundant in the 
our model, which leads to even stronger constraints. In other words, the application 
of dark Higgs results in ref. [76] here leads to the conservative constraints. 

• Cosmic Microwave Background (CMB): the CMB was formed when the photon 
last scattering occurred at ~ 10 12 s after the Big Bang. For most of the parameter 
space the duration of CMB formation 10 12 s) was much longer than the /12 lifetime. 
Thus annihilation of the VDMs into a pair of ft2 [77], which further decayed into ener¬ 
getic electrons and positrons, distorted the CMB spectrum in Cosmic Dark Ages [78- 
80]. Moreover, such a CMB constraint was further strengthened by the Yukawa poten¬ 
tial in eq. (4.3) via the Sommerfeld enhancement [81, 82]. For the s-wave VDM anni- 
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hilations, this can lead to a nonperturbative correction to the tree-level cross section, 
av = S x (crv) o, in which (cru)o denotes the tree-level perturbative cross section for 
XX /i 2^2 and S is the s -wave Sommerfeld enhancement factor given by [29, 83-85] 


S=* - 3inh(2TO<:> , ( 6 . 1 ) 

a cosh(27 vac) — cos(27ry c — (ac) 2 ) 

with a = v/(2ax) and c = 6axmx/(7r 2 rrih 2 ). Since the velocity of the VDM was very 
small during the photon last scattering, we can use the value of S saturated in the 
vanishing velocity limit. Due to the large mass hierarchy between the VDM X and 
the mediator / 12 , the CMB upper limit in figure 8 of ref. [ 86 ] for the one-step cascade 
with the e + e - final state can be applied for the VDM annihilation cross-section. 

• AMS-02: the local annihilations of VDMs into h<i pairs decaying to e + e - in the 
final state can lead to an excess of positron flux in cosmic rays [87-89]. Therefore, 
the absence of such an excess would give rise to a strong upper bound on the VDM 
annihilation cross section. Currently, the most precise measurements of the positron 
flux [90] and positron fraction [91] come from the AMS-02 Collaboration. By taking 
into the account the Sommerfeld enhancement factor in eq. (6.1) with typical VDM 
velocity vx ~ 10 -3 in our Galaxy, we can take the AMS -02 positron flux constraints 
from ref. [ 86 ] for one-step cascading VDM annihilations. Note that the AMS-02 
results are reliable only down to the DM mass ~ lOGeV, since the positron flux 
spectrum lower than 10 GeV would be affected significantly by the solar modulation 
so that the constraints in this range would be uncertain. 


• Dwarf limits from Fermi: the VDM annihilations in the dwarf spheroidal galaxies 

provide bright 7 -ray sources in the Milky Way, and are thus expected to be probed 
and constrained by the Fermi Gamma-Ray Space Telescope [92]. In the present 
model with > 2 m e , most 7 -rays are generated by the final-state radiation from 
the mediator decay e + e _ 7 , which follows the VDM annihilation XX ^ 2 ^ 2 • 

However, due to the suppression factor from radiative corrections compared with the 
dominant decay channel h<i —» e + e - , the constraints from Fermi shown in ref. [ 86 ] 
are much weaker than the corresponding ones from CMB and AMS-02. Therefore, 
we do not show dwarf limits from Fermi in our following numerical results. 

Region (ii): m^ 2 < 2 m e . We now turn to the indirect search constraints for the 
VDM with the mediator mass m ^ 2 < 2m e , in which decays dominantly in the diphoton 
channel, and the lifetime is typically longer than 10 12 s. As mentioned before, for such a 
light / 12 , the BBN constraints can be evaded as shown in ref. [76]. 

• Dwarf limits from Fermi: since /12 —>* 77 is the dominant /12 decay we expect that 
there should be strong constraints from measurements of 7 -rays by Fermi Gamma- 
Ray Space Telescope [92]. However, note that the signal region for each dwarf is 
defined as the one within an angular radius of 0.5°. For the 15 dwarfs used in the 
Fermi-LAT analysis, their distances from the Earth range from 32 kpc to 233 kpc. 
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Thus, due to the fact that propagates at the speed of light without any scatterings 
in the range of a dwarf, the h<i will spend, at most, the time of 0 (lO n s) traveling 
inside the signal region from the center of the dwarf. In other words, it is too short in 
time for to decay inside the signal region. As a result, the Fermi-LAT constraints 
in ref. [92] cannot be adopted directly in our case. 


• CMB: for r ^ 2 > 10 12 s, h<i would have a large abundance at the time of recombi¬ 
nation. Also, the high-energy photons from /12 decays would ionize and heat neu¬ 
tral hydrogen after recombination, and hence distort the CMB anisotropy spectrum. 
Consequently, recent measurements of the CMB by Planck [77] can provide strong 
constraints on properties [93, 94]. We adopt the recent lower bound on the decay¬ 
ing DM lifetime r° for the diphoton final state shown in figure 7 of ref. [95] to obtain 
the following constraint for the VDM model: 


Th 2 > r° x 


0.12 ’ 


(6.2) 


where D/^/z 2 is the current /12 relic density generated via the freeze-in mechanism 
if /12 were present today without decays. In fact, h<i might decay well before. The 
constraint is actually for abundance at the epoch of recombination, not today. 
We only use the present DM relic density as a reference to quantify the /12 density 
fraction at the recombination period. Moreover, the expression on the right-hand side 
in eq. ( 6 . 2 ) is just an approximation and the true formula should be D/ l 2 /i 2 /(D/ l 2 /i 2 + 
£lxh 2 ). However, the /12 density is always smaller than that of VDM due to the 
assumed mass hierarchy, so that D^/i 2 in the denominator can be neglected. Note 
that the exclusion limit in ref. [95] extends to the DM mass of lOkeV, so we ignore 
the CMB constraints below this VDM mass in our numerical calculations. 


• Diffuse 7 /X-ray bounds: when the lifetime of h<i is larger than the present age 
of the Universe tjj = 4.3 x 10 17 s, the particle contributes to the present DM relic 
density even though it is not absolutely stable. The only decay channel h<i —> 77 
could be constrained by the accurate measurement of the diffuse 7 /X-ray background. 
Following refs. [44, 96-98], we adopt the conservative lower limit on the /12 lifetime as 

Th 2 > 10 28 sec x , (6.3) 

where is the relic abundance of generated via the freeze-in mechanism. 


7 Numerical results 

Having discussed all of the VDM signals and constraints, we can put everything together 
to see if we can find a region in the parameter space where large DM self-scatterings for the 
scale of dwarf galaxies can be compatible with the VDM relic density and all of indirect 
search constraints. Note that there are four free parameters in our original VDM model, 
so that if the Higgs portal coupling k is fixed as shown in figure 4 by the requirement that 
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m h = 100 MeV 



1 2 
Log 10 (m x /GeV) 



Figure 7. Region (i) in the parameter space with m^ 2 = 100 MeV (left panel) and rrih 2 = 2 MeV 
(right panel) for which the VDM self scattering might be sufficiently strong to solve the small 
scale structure anomalies. The numbers next to “dw” represent crr/mx in units of cm 2 /g at the 
dwarf scale. The blue (cyan) region shows the space with the VDM self-interaction cross section 
of 1 cm 2 /g < (Tt/ttix < 10 cm 2 /g (0.1 cm 2 /g < (Jr/mx < 1 cm 2 /g) at the dwarf scale, while 
the red region shows bounds on VDM self scatterings at the cluster scale. The purple, green, and 
orange colors show regions excluded by the DM indirect searches from BBN, CMB, and AMS-02, 
respectively. The region below the curve named “Thermalization” shows the parameter space in 
which the VDM is generated via freeze-in, and the gray region above the curve corresponds to the 
one with dark sector thermalized in itself. In both plots the Higgs portal coupling k is fixed by the 
VDM relic density as a function of mx in the freeze-in region. 


the VDM relic density constitute all of the DM in the Universe, we can plot the parameter 
space in the mx-&x plane with fixed values of m^ 2 . The final results for some typical 
values of rrih 2 in the Regions (i) and (ii) are presented in figure 7 and figure 8, respectively. 

Figure 7 shows constraints on the parameter space for two extremal values of rrih 2 
in the Region (i). Note that only the region below the curve named “Thermalization” 
represents the freeze-in region, so that above the curve the parameter space corresponds 
to the case with dark sector thermalized by itself. It is seen that the DM indirect search 
constraints from BBN, CMB, and AMS-02 are so strong that the whole freeze-in region 
is excluded, no matter how precisely we tune the model parameters. If we focus on the 
region that solves the small-scale structure problems, the left panel shows that the only 
parameters that exactly sit at the resonances can enter the freeze-in region for a relatively 
heavy h<i (rrih 2 = 100GeV). In this case, the most stringent constraints come from the 
CMB and AMS-02, which can be understood that the same resonances that give rise to 
the appropriate VDM self-interactions would also induce large Sommerfeld enhancements 
in the VDM annihilations. However, as the /i 2 mass decreases, more and more signal 
regions at the dwarf scale are shifted to the classical and Born regions. In the extremal 
case with mh 2 — 2 MeV which is chosen to avoid the e + e - threshold effect around 1 MeV, 
the velocity dependence of the VDM self-scatterings becomes manifest, as the constrained 
region from clusters separates from the signal region at the dwarf scale. However, the whole 
signal region below the thermalization curve is excluded by the combination of the CMB, 
AMS-02 and BBN constraints. 
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m h = 0.1 MeV 



1 2 
Log 10 (m x /GeV) 



Figure 8. Same as in figure 7 but for Region (ii) with rrih 2 = 10 “ 1 MeV (left panel) and rrih 2 = 
10 -3 MeV (right panel). When the lifetime of h 2 is longer than the current age of the Universe 
tjj = 4.3 x 10 17 s (denoted by the dashed black curve), we also include the constraints from the 
diffuse 7 /X-ray observations. 


The situation changes a lot for the Region (ii) as shown in figure 8 , since the indirect 
detection constraints are all imposed on the decay process 77, rather than VDM 

annihilations. Here we only consider the freeze-in region below the thermalization curve. 
In both panels, the signal regions for the dwarf galaxy scale are all in the Born and classical 
regions, in which the part with a small hidden gauge coupling ax and a light VDM corre¬ 
sponds to the Born region while the band with large values of ax and mj to the classical 
region. The discontinuities in both plots represent the mismatch of the analytical formulae 
around the boundary of these two regions. By detailed calculations, it is found that all of 
the signal regions for rrih 2 > lO -2 MeV are constrained tightly by observations of CMB 
and diffuse 7 /X-rays, as illustrated by the left panel of figure 8 . Only when the mass 
is reduced to O(keV) a small parameter window opens, in which the VDM mass is around 
O(GeV) and ax is in the range 10 -9 ^ 10 -6 , as is seen clearly in the right panel of figure 8 . 

8 Conclusion 

We have investigated a simple VDM model in which the dark sector consists of only two par¬ 
ticles, the VDM X and the dark Higgs / 12 , and it couples to the SM sector weakly through 
the Higgs portal. We are particularly interested in the region of the parameter space with 
the VDM mass of 0(GeV ~ TeV) and the /12 mass of m ^ 2 < 100 MeV, where the dark 
Higgs /i 2 plays a role of a light mediator so that non-perturbative effects can generate VDM 
self-interactions with the appropriate magnitude to solve the small-scale structure problems 
at the dwarf galaxy scale. Due to the velocity dependence, such VDM self-scatterings can 
avoid the constraints at the galaxy cluster scale. In our work, we consider the scenario in 
which the VDM relic density is produced by the freeze-in mechanism. Especially, we pay 
attention to consequences of the EW phase transition, since they affect the dominant DM 
production channels quite dramatically. It turns out that the Higgs portal coupling k is al- 
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ways predicted to be of 0(1O -11 ). With such a tiny portal coupling, the DM direct detection 
limits cannot constrain the model much, as it is illustrated by the latest complete LUX data. 

However, the indirect detections can place strong constraints on the parameter space of 
interest. Specifically, when 2 m e < m^ 2 < 100 MeV, decays dominantly via h<i e + e~ 
so that the limits from BBN, CMB, and AMS-02 exclude all of the parameter space. In 
particular, the CMB and AMS-02 bounds are strengthened since the main VDM annihila¬ 
tion XX —>► /i 2^2 is enhanced by the non-perturbative Sommerfeld effects. On the other 
hand, if < 2 m e , the relevant indirect search constraints come from observations of 
CMB and diffuse 7 /X-rays. As a result, only when the mass of is equal to or smaller 
than keV scale one can find a parameter window in which the model can lead to the right 
DM relic abundance and appropriate DM self-interactions, while it does not conflict with 
other indirect detection observations. 
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A Relevant cross sections and decay rates for the Higgs portal 


In this appendix, we collect the formulae for the relevant SM particle annihilation cross 
sections as well as SM-like Higgs hi decay rate, which are involved in the calculation of dark 
matter relic density in the Universe via the freeze-in mechanism. Since the SM EW phase 
transition has a substantial impact on the VDM production, we consider the annihilation 
and decay channels in broken and symmetric phases, respectively. 


EW symmetry-broken phase 
Quark annihilation to VDMs: 


a(qq -> XX) = 


K 2 m 2 

(s — 4m 2 )(s — 4m 2 x ) 

s 2 — 4m x s + 12m^ 

192vr \ 

/ s 2 

(s - m 2 2 ) 2 [(s - m 2 hi ) 2 + r 2 h m 2 hi \ 


(A.l) 


where we use the SM Higgs boson width r h = 4.15 MeV [99, 100] to regulate the SM-like 
Higgs mass pole singularity. 


Lepton annihilation to VDMs: 


cr(ll -> XX) = 


K 2 m 2 / (s — 4m 2 ) (s — 4m x ) 


s 2 — 4m 2 x s + 12mt 


X 


(j{W + W~ -> XX) = 


647 TV S 2 ( S “ m D 2 [( S - m hi) 2 + r hi m h 1 ] ' 

hi 2 I (s — 4 rriyy)(s 4m x ) s 2 — 4m x s + 12m^ 


(A. 2 ) 


W bosons annihilation to VDMs: 

.2 


2887T 


0 


m 




>[( S 


m 


h\ > 


+n 




x- 


s 2 — 4m^ v s + 12 m 


w 


s - 4 


(A.3) 
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Z bosons annihilation to VDMs: 


(j{ZZ -> XX) = 


ft 2 / (s — 4m 2 z )(s — 4m x ) 


288?r 


s 2 — 4m 2 x s + 12m 4 


•x 


( s - - m l) 2 + r l m l ,] 




/u ##< 7uJ 


s 2 — 4 m 2 z s + 12m z 
s — 4 m 2 z 


(A.4) 


SM-like Higgs bosons hi annihilaiton to VDMs: 


Tv 2 

_ A 

7/T 

1 

7/T 

1 

s 2 — 4m 2 x s + 12m^ 

“ 32t r \ 

1 s 2 

(s - m 2 2 ) 2 [(s - m 2 hi ) 2 + T 2 h m 2 hi ] 


2 \2 


( s + 

s - 4m 2 hi 


(A.5) 


where we have only kept the leading-order terms in the double expansion of ft and ra| 2 /ra^. 

SM-like Higgs boson hi decay to VDMs: when the VDM mass is smaller than a 
half of the SM-like Higgs mass, the VDM can also be produced by the decay of hi, with 
the decay rate as follows: 


r(hi -> xx) = 


"h! 


2k 2 ^ 2 c 2 


X m h 1 
,2 ™4 


K, - 2 m 2 x) 2 
4m x 


2 + 


a 


( m h - m L)V 


— 4m^- 


(A. 6 ) 


( m fe t ~ 2m^) 2 

4m^ 


2 + 


where we have used the definition of ft in eq. (2.7) and the approximation that cq « 1. 

SM-like Higgs boson hi decay into h^’s: 


r(hi —>• /12/12) 


.2„,2 


K Vrr 


32nm hl ’ 


(A.7) 


where we only keep the leading order in the expansion of ft and m/ l 2 /m/ ll because m ^ 2 <C 
m hl . 


EW symmetric phase 

SM Higgs doublet H annihilation to VDMs: 


a(XX -y HH ] ) 


k 2 I (s — 4 rrifj)(s — 4m|-) s 2 — 4777,^ s + 12 tt7 
727 rV s 2 (s — m s ) 2 (s — 4m 2 x ) ’ 


(A. 8 ) 


where tuh and m s are the masses of the SM-Higgs doublet and the dark Higgs (j)s defined 
in eq. (2.4). 

Note that in our derivation of the Boltzmann equation in eq. (3.1) after and before the 
EW phase transition, we have used the so-called reaction density 7 % for various channels. 
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For the annihilation channels, the reaction density is defined as [101] 


7 (a&-» 12 ) = J dp a dp b dp 1 dp 2 fa <i fb q (‘2Tr) 4 5 4 (p a + p b - p 1 - p 2 )\M{ab -> 1 2) | 2 (A.9) 

roc 
>1 . 


32vr 4 


r °° 


V~s 


12)^1 


where a, b (1,2) represent the incoming (outgoing) particles with g a ^ as their respective 
degrees of freedom, and f^ q ~ e Ei/T 

is the Maxwell-Boltzmann distribution. Here dp = 
d 3 p/(27r) 3 (2E), \M\ 2 are the amplitude squared summed over quantum numbers of the 
initial and final states without averaging, and s m i n = max[(m a + m^) 2 , (mi + m 2 ) 2 ]. 

On the other hand, the reaction density for the decay channel a -A 12 has the following 
definition [38, 41]: 

7 D (o ->• 12) = J dpadp 1 p2(27r) A 8 4 (pa - Pi -p 2 )fa q \M(a ->• 12) I 2 

= J^r(a 1 2)Ti£ 1 (^) , (A.10) 

where is the zero-temperature decay rate. 
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